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Abstract 


Multibody mechanical systems find a tremendous number of applications in 
various fields like I. C. engines, automobile transmission, robotics, structural 
supports etc. Analysis of such systems leads to massive algebraic manipula- 
tions in constructing equations of motion. Since these equations are highly 
nonlinear, the prospect of obtaining closed-form solutions is remote, except 
in very simple cases. The objective of computer-aided analysis is to formulate 
and solve such problems using digital computer. Several computer programs 
for analysis of mechanical systems are developed since 1970’s. 

In the present work, the static equilibrium analysis of compliant mechani- 
cal systems under the action of external loading is attempted. Statically de- 
terminate as well as statically indeterminate systems are analysed. Relative 
coordinates are used as generalised coordinates. Using relative coordinates 
not only requires less number of generalised coordinates, but it also makes 
analysis computationally efficient. The system is modeled using homoge- 
neous transformation matrices and loop closure equations. For the solution 
of the equilibrium problem, potential function (<I>) is used. Minimisation of 
<1>, subject to rigid constraints gives equilibrium configuration. The approach 
is tested on some planar and spatial mechanical systems. 
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Chapter 1 
Introduction 


A mechanical system is defined as a collection of bodies (or links) in which 
some or all of the bodies can move relative to one another (Fig. 1.1). Broadly 
they are classified into two types, mechanisms with degrees of freedom greater 
than zero and structures, with zero degree of freedom. 

Mechanical systems find wide range of applications in various fields. A 
few examples are automobile transmission, robotics, structural supports, I. 
C. engines, suspension systems, aviation etc. 

1.1 Simulation of Mechanical Systems 

The process which allows the engineer to study the response of an already 
existing system to a known excitation is called analysis. This requires a 
complete knowledge of the physical characteristics of the mechanical system, 
such as material composition, shape and arrangement of parts. Due to the 
nonlinear nature of analysis problems, mechanical designer has traditionally 
resorted to graphical techniques and physical models for system analysis. As 
might be expected, such methods are limited in generality and rely heavily 
on the designer’s intuition. With the advent of digital computers, engineers 
started using the computer and available numerical methods to solve highly 
nonlinear analysis problems. 

Factory automation is one of the major objectives of modern industry. 
Various parts of the product are designed in the Computer-Aided Engi- 
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neering branch. Computerised product design requires such capabilities as 
Computer-aided analysis (CAA), Computer-aided drafting, design sensitivity 
analysis and optimisation. The Computer-aided analysis capability serves as 
part of the design process. It allows the engineer to simulate the behaviour of 
a product and modify its design prior to actual production. Thus Computer- 
aided analysis can be used as a partial substitute for laboratory or field tests 
in order to reduce the cost. 



Figure 1.1: Mechanical systems with compliant members 


3 



1.1.1 Advantages of computer simulation 

Computer simulates the behaviour of various types of multibody systems 
in great detail during design stages, from first design concepts to the final 
prototypes. At any stage of design process, CAA provides an important 
amount of data to study the influence of the different design parameters, 
even before the construction of any physical device or prototype. These 
tools allow the designer to carry out a large number of simulations quickly 
and economically. Thus mechanical system analysis can be made more useful 
to the designer by having the computer carry out repetitive calculations. 
The engineer can then concentrate on more creative concepts/aspects of the 
design process. 

1.1.2 Requirements 

A general purpose computer program for the analysis of mechanical systems 
must perform following basic functions. 

• Accepting data from the user. 

• Generating the governing equations of motion. 

• Solving the equations automatically. 

• Communicating the result to the user. 

1.2 Static Equilibrium Analysis Problem 

Under the name “ position problem ” fall several different, although closely 
related, problems. 

1.2.1 Initial position problem 

The initial position problem is concerned with the determination of the po- 
sitions of all links of the mechanism, knowing the position of the input links 
(as many as there are degrees of freedom in the mechanism). Solution of 
the position problem is important because often it is a preliminary condition 
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for kinematic and dynamic analysis. In this thesis, this problem is referred 
to as assembly phase, where from approximate input given by user, correct 
configuration of the system is calculated by using optimisation technique. 

1.2.2 The problem of finite displacement 

This problem is closely related to the former problem and it lies in finding 
the displacement of the different links for certain displacements of the input 
links. The problem can be approached by the same nonlinear techniques as 
the initial position problem. 

1.2.3 Static equilibrium position problem 

Finally, this problem determines, as its name indicates, the equilibrium posi- 
tion of the mechanical system subjected to known external forces, including 
the ability to store potential elastic energy in its compliant members. Fig. 

1.2 shows initial configuration of spring restrained slider crank mechanism 
in thin lines. After applying force F, equilibrium configuration of system is 
shown in thick lines. 

Unlike velocity and acceleration analysis, above mentioned position prob- 
lem is strongly nonlinear and for this reason, has always presented much 
greater difficulties for its solution. Kinematic and dynamic analysis of com- 
plex mechanical systems is often initiated from a position of static equilib- 
rium. Assigning correct values to the coordinates that describe a state of the 
static equilibrium, can be a complicated task for a large-scale inter-connected 
systems of bodies. 

1.3 Literature Review 

Some methodologies suggested to solve static equilibrium problems are dis- 
cussed in this section. 
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Figure 1.2: Static equilibrium for simple mechanical system 

1.3.1 Minimum P.E. and structural strain method 

The work of Livermore [1] is considered classic in this field. He describes 
two methods developed by himself and by Gupta [2], called, respectively, 
“Minimum Potential Energy Configuration Method’ and “Structural Strain 
Method’ . The first process is slow and laborious, especially when the system 
has several degrees of freedom. The second method begins with a known 
equilibrium position at which the mechanism has no external forces applied. 
Then, an incremental process begins, in which external forces are applied 
little by little and the geometry of the mechanism changes as the forces are 
applied. This method is based on the analysis of displacement and veloci- 
ties. Livermore points out that the calculation time according to the second 
method is usually between 1/2 and 3/4 the time required by the first. But it 
requires more complex code and considerably greater storage capacity. Also 
the method may converge to unstable equilibrium, so the stability of solutions 
must be checked. 
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1.3.2 Gas method 


Most attractive features of this method [3, 4] are conceptual simplicity and 
the ease with which it can be programmed. It offers a coherent unified way to 
approach kinematic, static equilibrium, kinetostatic and dynamic problems. 
It uses basic coordinates. These are the coordinates of some previously chosen 
points of the mechanical system. The problem is solved like a problem of 
mathematical programming, and a minimisation method is proposed. 

1.3.3 Equilibrium position from equations of motion 

Alternative method where static equilibrium position is achieved by setting 
velocity and acceleration vectors to zero in equations of motion is suggested 
in [5, 6]. But resultant algebraic equations are highly nonlinear and initial 
estimates of some of the unknowns, such as Lagrangian Multipliers are not 
readily available. Also poor estimate of Lagrangian Multipliers may lead to 
badly conditioned matrix and divergence of the algorithm. 

1.3.4 Force equilibrium approach 

The basic methodology of this approach is described in [7]. In it’s detailed 
implementation [8], the work is divided into two parts. In the first part, on 
the basis of approximate information about the system, provided by the user, 
the initial correct configuration of the system is obtained. In the second part, 
the equilibrium analysis of the system under the given loading is carried out 
to get the final equilibrium configuration. 

The basic methodology used is to develop the governing constraint equa- 
tions for rigid and deformable members. These equations are then used in 
the force equilibrium equation to get the deformed configuration of the sys- 
tem at equilibrium. In this approach, generalised coordinates used are the 
relative coordinates of every link with respect to the previous link in the 
loop. Deformable members are modeled as linear and torsion springs. 

Although above suggested method is computationally fast, accurate and 
efficient, it has serious limitations, since it only deals with statically deter- 
minate systems, which severely restricts its application area. In this thesis, 
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the above approach of modeling the system using relative coordinates is used 
to develop a general formulation which can analyse statically indeterminate 
systems also. 

1.3.5 Potential function approach 

In this approach [9], modeling of mechanical system is based on construction 
of a simple F.E.M. with truss elements. For solution of the problem, error 
function (<&) is the potential function of the mechanical system. This func- 
tion is the sum of two terms, the elastic potential of the deformable elements 
and the potential of the set of conservative external forces transformed to 
the nodes of the model. The error function is subjected to constraints that 
lengths of the rigid elements remain constant. Solution is sought through 
various optimisation schemes and the methodology covers all position prob- 
lems. 


1.4 Scope of the Present Work 

The thesis deals with the static equilibrium of user defined planar and spatial 
mechanical systems with compliant members. The mechanical system can 
be statically determinate or statically indeterminate. 

After the user provides approximate input of the system, correct initial 
configuration is obtained by imposing constraint equations. Then equilibrium 
analysis of the system for given loading is carried out using minimisation of 
potential energy. 

1.5 Assumptions 

• Compliant members are considered as linear elastic elements, i.e. there 
is a linear relationship between the forces applied and the deformation 
of these members, irrespective of magnitude of forces and deformation. 

• Only binary flexible links are considered. Rigid links can have any 
number of connections. 
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• Friction due to relative motion of links at joints is neglected. Also the 
damping present at joints and in the links is neglected. 

• Joint clearances are neglected. 

• The loading applied on the system is time independent and stationary, 
i.e. the absolute position of the forces remain the same during the 
deformation. 

1.6 Organization of the Thesis 

• Chapter 2 : It introduces general formulation of the problem. 

• Chapter 3 : Modeling of systems with various joints is described in 
detail in this chapter. 

• Chapter 4 : It describes assembly and static equilibrium analysis pro- 
cedure. 

• Chapter 5 : In this chapter, some representative results of assembly and 
static equilibrium analysis of planar and spatial mechanical systems are 
reported to establish the capability of the method. 

• Chapter 6 : Conclusions and future scope is presented in this chapter. 



Chapter 2 
Formulation 


This chapter introduces various types of generalised coordinates with their 
relative comparison. In later part of the chapter, general formulation of static 
equilibrium problem is explained. The methodology is based on [7, 8, 9]. 

2.1 Generalised Coordinates 

Any set of variables that uniquely specifies the position and orientation of 
all bodies in the system, that is, the configuration of the system, can be 
taken as the set of generalised coordinates. The computational efficiency of a 
general-purpose program depends upon several factors, two of which are the 
choice of coordinates and the method of numerical solution. The choice of 
coordinates directly influences both the number of equations of motion and 
their order of nonlinearity. 

2.2 Selection of Generalised Coordinates 

In order to show how different sets of coordinates can lead to different for- 
mulations describing the same system, a simple example is given here [5]. In 
this example, the four-bar linkage is considered for kinematic analysis. 

The first formulation shown here considers only one coordinate to describe 
the configuration of the system. This is referred to as the minimal coordinate 
of the system. As shown in the Fig. 2.1(a), the angle (j), describing the 
orientation of the crank with respect to the X axis, can be selected as the 
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Y 




(b) Relative coordinate 


Y 



Figure 2.1: Selection of generalised coordinates 
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generalised coordinate. For any given configuration, i.e. known 0, any other 
information on the position of any point in the system can be calculated l . 

The second method considers three coordinates 0i, 0 2, 0 3. Fig. 2.1(b) 
shows the arrangement. For this mechanism, the selected coordinates de- 
fine the orientation of each moving body with respect to previous body. 
Therefore, these coordinates are referred to as relative coordinates. Since the 
four-bar linkage has only one degree of freedom, the three coordinates are not 
independent. Two loop equations relating these coordinates can be written 
as follows. 


T COS (0i) + d COS (0i + 02) + S COS (0i + 02 + 03) — l — 0 

r sin (01 ) + d sin (0! + 0 2 ) -I- s sin (0i + 0 2 + 0 3 ) - v = 0 

The third formulation uses three Cartesian coordinates per link, the x and y 
coordinates of the center point of each link and the angle of the link which 
is measured with respect to the X axis. The arrangement is as shown in the 
Fig. 2.1(c). So there are total nine coordinates dependent upon each other 
through eight equations. 

The crude but general comparison between these three sets of coordinates, 

with regard to several crucial and important aspects, is summarized in table 

2.1 2 . A general conclusion that can be made from this table is that the 

smaller the number of coordinates and equations, the higher the order of 

nonlinearity and complexity of the governing equations of motion, and vice 

versa. From above comparison it is clear that a set of relative coordinates 

falls in the middle of the “comparison scale”. Therefore the selection of a 

set of relative coordinates is a good compromise between other two extremes, 

and the same is used in our formulation. It is also useful in getting derivatives 

of the governing equations in a simplified way. These derivatives arc required 

in analysis. 

hipto the branch 
included from reference [5] 
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Minimal 

coordinates 

Relative 

coordinates 

Cartesian 

coordinates 

Number of coordinates 

Minimum 

1 

Moderate | 

1 

Large 

Number of second-order 
differential equations 

Minimum 

Moderate 

! Large 

| 

Number of algebraic 
constraint equations 

None 

Moderate 

Large 

Order of nonlinearity 

High 

Moderate 

Low 

Derivation of the 
equations of motion 

Hard 

^ Moderately 
hard 

Simple 

Computational 

Efficiency 

Efficient 

Efficient 

Not as efficient 

Development of a 
general-purpose 
computer program 

Difficult 

Relatively 

difficult 

Easy 


Table 2.1: Comparison of generalised coordinates 


2.3 Constraint Equations 

In mechanical system, a kinematic pair imposes certain conditions on the rel- 
ative motion between the two bodies it comprises of. When these conditions 
are expressed in analytical form, they are called equations of constraint. A 
constraint can be seen as any condition that reduces the number of degrees 
of freedom in a system. The algebraic equations of constraint must be de- 
rived such that they are equivalent to physical joints. This equivalence is 
important and constraint equations must imply the geometry of the joint. 
In this step, the constraints (if any) among the generalised coordinates are 
developed in the form as shown below. 
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feR m 


( 2 . 1 ) 


f(X) = 0 

Where, m = number of rigid constraints 
i.e. 

/t(0t) 02) 03) ~ 0 

02(01) 02) 03) •••0n) = 0 


/m(0 1) 02) 03) 0n ) — 0 

Where, X = [0i, $ 2 , 03) is the set of generalised coordinates. 

2.4 Static Determinacy 

As the mechanical system under consideration is constrained and its mobility 
is owing to the deformable members only, we have 

m + p > n 

Here p is number of deformable coordinates and m is number of rigid con- 
straints.. The cases for which the equality holds are statically determinate 
cases and those for which the strict inequality holds are the statically indeter- 
minate cases. Fig. 2.2(a) shows a statically determinate system comprising 
of three rigid and two deformable members. In Fig. 2.2(b), one additional 
torsional spring at base A, make the configuration statically indeterminate. 
In a planar system, every link has three degrees of freedom. In Fig. 2.2(a), 
there are three rigid links, of which one is ground link. Hence, for remaining 
two links we have n = 6. Also p = 2, due to two deformable members (links 
2 and 3). For constraints, we have to consider various joints in the system. 
Every revolute joint connecting two rigid links gives two constraints. Due to 
two revolute joints (Ri and R 2 ), rri is 4. Since m + p = n, this is statically 
determinate system. 

In Fig. 2.2(b), due to additional torsion spring, p becomes three. In this 
case {rn + p > n), and system is statically indeterminate. Similarly, the 
mechanical system shown in Fig. 2.2(c) is also statically indeterminate. 
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Figure 2.2: Statically determinate and indeterminate cases 

2.5 Generalised Force 

For the static equilibrium, forces are considered to be acting at the fixed 
locations and along the fixed orientations, with respect to the global reference 
frame, throughout the deformation of the system. The loading constitutes of 
the external stationary forces and moments acting on the system. Weights 
of various links can also be considered as a part of loading. Input force is 
transformed to generalised force by calculating Jacobian. 
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2.6 Static Equilibrium Analysis 

Here static equilibrium configuration of the mechanical system is determined 
by evaluating the position for which the potential energy (P.E.) function is 
at its minimum. The factors constituting the system P.E. are, gravitational 
field, externally applied forces and compliant members. 

We have used $ as the P.E. function of the mechanical system. This 
function is the sum of two terms, 

• the elastic potential of the deformable elements and 

• the potential of the set of conservative external forces transformed to 
the nodes of the system. 


*({*}) = \ E M4(W) - A) 2 + i £ M - {F } T {AX} (2.2) 

z i=l A j = 1 

Here, 

* a is the number of linear springs 

* d{ is the length of the elastic element i in a given iteration 

* Di is the length without deformation 

* ki and k t j are stiffnesses of linear and torsional springs respectively 

* {F} : generalised force 

* b is the number of torsion springs 

* A 9 is the angular deformation of torsion spring 

* {X} is the generalised coordinate vector 

However, since the mechanical systems can include both rigid and elastic 
elements, it is necessary that the lengths of the rigid elements remain con- 
stant. In other words certain constraints must be satisfied, as discussed in 
section 2.3. The number of constraints is decided by factors such as whether 
the system is planar or spatial, number of loops, and type of joint at the 
end of the loop. Minimisation of potential function in Eq. 2.2, subject to 
constraints in .Eq. 2.1, gives the final equilibrium configuration. 
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Chapter 3 

System Modeling 


A closed loop mechanical system may contain one or more loops. Present 
work can analyse single as well as multiloop systems. Every loop starts from 
the ground connection of the link at which global reference frame is fixed and 
ends at the another ground connection. Position of other ground connections 
in global reference frame should be fixed. 

3.1 System Input 

The set of input required for describing the mechanical system is : 

• Nature of the system ( planar or spatial) 

• Total number of links in the system 

• Number and type ( linear or torsion ) of deformable members 

• Connectivity of links in a loop 

• Link dimensions 

• Ground connection coordinates with respect to global reference frame 

• External loading ( forces , moments ) with their location and orientation 

• Type of joint at ground connection 
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• Approximate values of the relative coordinates of the links with respect 
to previous link along the loop . 

The order in which this input is to be provided by the user, is given in the 
appendix (A), for a sample problem. It is important here, to note that the 
user has the freedom to enter approximate values of generalised coordinates 
(g.c.) only. The other input parameters such as link lengths and ground 
connection coordinates must be provided as accurately as possible. 

3.2 Link and Joint Modeling 

The mechanical system is modeled using concatenated homogeneous trans- 
formation matrices as discussed in Craig [11]. These matrices represent geo- 
metric relation of every joint in the system. 

For every link, local coordinate system is attached. Fig. 3.1 shows two 
links with their coordinate systems. They are connected by revolute joint. 
The x-axis of every local coordinate system is along the direction of connec- 
tion between the two links. The links revolve about the 2 -axis. 

3.2.1 Revolute joint 

It offers one degree of freedom. For a planar case shown in Fig. 3.1, link i 
and revolute joint ij are described by using translation of L — [l x ,l y , l z , 1], 
from local reference frame of link i to that of link j. Subsequently, there is 
rotation about 2 axis through angle 0. Therefore, the relation of the j-th 
link with respect i-th link is 

T rev = T[L ) x A(0) 


■ 1 

0 

0 

l x 1 


COS 0 

— sin 0 

0 

0 ‘ 

0 

1 

0 

h 


sin 0 

COS 0 

0 

0 

0 

0 

1 

lz 


, 0 

0 

1 

0 

. 0 

0 

0 

1 . 


0 

0 

0 

1 . 
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Figure 3.1: Revolute joint ij 

cos <fi —sin <f> 0 l x 
sin <f> cos <j> 0 l y 

0 0 1 l z 

0 0 0 1. 

In this case l x = l y — 0. Thus, T rev (L, <f>) is the homogeneous transformation, 
describing the revolute joint. L is the translational distance and (f> is the 
amount of rotation with respect to previous link. With this transformation 
matrix we can express the coordinates on the j- th link with respect to the 
i-th link. 

In case of spatial mechanical systems with revolute joint, the procedure is 
as follows [10]. Suppose (a, b, c) are direction cosines of axis and (xi,yi,zi) 
is one point on it. First we reposition the axis so that it passes through the 
origin. Transformation matrix for it is given by, 

'10 0 —x\ 

T= 0 10 -^ 

0 0 1 -z x 
. 0 0 0 1 

Then, we make rotation axis parallel to Z axis by first rotating it about X 
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axis and then about Y axis. 


Here, 


R x {a) 


Rv(0) 


0 0 


0 c/d 

-b/d 0 

0 b/d 

c/d 0 

0 0 

0 1 

' d 0 

-a 0 ' 

0 1 

0 0 

a 0 

d 0 

0 0 

0 1 




and 


sin(o;) = b/d cos(a) = c/d 


Now, we have aligned the rotation axis with the positive Z axis. The specific 
rotation angle 9 can now be applied as a rotation about the Z axis. To 
complete the required rotation about the given axis, we need to transform 
the rotation axis back to its original position. This is done by applying 
inverse of above transformations. The transformation matrix of rotation 
for revolute joint then can be expressed as the composition of the several 
individual transformations : 


T rev = T~ l * R~ l {a 0 * R; l (f3) * R z (9 ) * Ry(/3) * R x {a) * T (3.1) 


3.2.2 Prismatic joint 


Here relative displacement between the two links connected is A L along z 


axis. The procedure is same as revolute joint, except that R z {9) in Eq. 3.1 
is replaced by T pr (AL ) as shown below. 


Tp r (AL) = 


'10 0 
0 1 0 
0 0 1 
0 0 0 


0 

0 

AL 

1 


(3.2) 
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3.2.3 Cylindrical joint 

It offers two degrees of freedom, translation as well as rotation along specified 
axis. Let A L be relative displacement between links and cp is relative angle 
between them. Then 


r C j,i(AL, <t>) = T(AL, 0) x R{ 0, (j>) 


Tcyl{AL, <p) — 


■ 1 

0 

0 

0 ■ 


" COS (p 

— sin 0 

0 

0 ' 

0 

1 

0 

0 


sin cf) 

COS (j) 

0 

0 

0 

0 

1 

AL 


0 

0 

1 

0 

. 0 

0 

0 

1 


0 

0 

0 

1 . 


(3.3) 


In case of spatial mechanical systems, matrix R( 0 , (p) is calculated similar to 
the procedure explained in revolute joint. 


3.2.4 Spherical joint 

It allows three degrees of freedom, i.e. rotation around all the three axis 
of one link relative to another. Here the Z-Y-X Euler angles are used to 
describe rotations. So, 


T sph (L, 4>z, <t> y , <j>x) = T{L, p z ) x R(4> y ) x R(<f> x ) (3.4) 


cos<fi 

— sin <j) z 

0 

L - 


COS 4>y 

0 

sin cpy 

0 ‘ 


' 1 

0 

0 

0 ' 

sin (f) z 

COS (j) z 

0 

0 


0 

1 

0 

0 


0 

COS (px 

“ sin cp x 

0 

0 

0 

1 

0 


— sin 4> y 

0 

COS (j) y 

0 


0 

sin <p x 

COS cpx 

0 

0 

0 

0 

1 . 


0 

0 

0 

1 J 


. 0 

0 

0 

1 . 


3.2.5 Screw joint 


It is similar to the cylindrical joint, except translation and rotation are related 
through pitch of screw. , The screw joint is represented by, 


T sc (/c, (j P) 


cos cp — sin <p 0 
sin <p cos (p 0 

0 0 1 

0 0 0 


0 

0 


L+g 


(3.5) 


Where, (p is the relative angle between the links, and k is lead of the screw. 
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Figure 3.2: Four bar mechanism 

3.3 Loop Closure Equations 

The governing equations are developed using loop closure equations. It con- 
sists of concatenation of homogeneous transformation matrices, along each 
loop in the system. A four-bar mechanism is used as an example to ex- 
plain the loop closure equation. It is shown in Fig. 3.2. Here, (X,Y,Z) is 
the global reference system. Links 1, 2 and 3 have local coordinate system- 
s (xi,yi,z l ),(x 2) y 2 ,z 2 ) and (x 3 ,y 3 ,z 3 ) respectively, fixed to them. Conse- 
quently, the (xi,yi,zi) coordinate system is defined relative to the global 
coordinate system by 

T t =T{ 0,&) 
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where is defined in the system 0. Furthermore, the (x 2 , ^ 2 ,^ 2 ) coordinate 
system is defined relative to the (aq ,yi,z x ) system by 

T 2 = T(Ly,fa) 

where L\ and 0 2 are defined in system 1. In a similar fashion, 

T 3 = T{L 2 , fa) 

where L 2 and 0 3 are defined in system 2. However, this does not close the loop 
for the mechanism. But, this matrix multiplication is sufficient for getting 
constraints, since ground connection is reached. The accurate location of end 
ground connection in this loop is already known. Let it be D = [x, y, z, 1] T . 
Then, constraint equations can be written as 

/ = T( 0, fa) * T(L U fa) * T(L 2 , fa) * {L 3 } - {D} = {0} (3.6) 

This gives two constraint equations (for x and y), in case of planar mechanical 
systems. Also fa, fa, fa being the generalised coordinates, these constraints 
are expressed as functions of them. 

Jacobian calculation, which needs differentiation of concatenated trans- 
formation matrices, is made easier by use of relative coordinates. It can be 
observed that, any one g.c. of the system described by loop closure equation, 
is wholly contained in only one of the matrices making up the product. Thus, 
while taking the derivatives of loop closure equations with respect to a g.c., 
say 4> 1 , only the components in the matrix having 0i are changed and the 
rest of the matrix product remains unchanged. 

Thus, we have, 

A- = Apr(0, 00 * T(L U fo) * T(L 2, fa) » {£ 3 }] 

<901 dfa 

= T' (0, fa) * T(L U fa) * T{L 2 , fa) * {L 3 } (3.7) 

Similarly, the derivatives with respect to all the g.c. can be calculated. 
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Chapter 4 

Static Equilibrium Analysis 


Analysis is done in two stages. In the first part, from the approximate input 
provided by the user, correct initial configuration is obtained. It satisfies all 
the constraints. This provision is useful in case of complicated planar and 
spatial systems, where giving accurate input system data is very difficult. In 
the next phase, static equilibrium analysis of the assembled configuration is 
performed and the final configuration is displayed in graphical format. 

The entire procedure is explained using six link planar system as shown 
in Fig. 4.1. It has rigid as well as deformable links. Whenever the system 
differs from this example (as in spatial cases), the difference and the relevant 
analysis is explained at suitable places. 

4.1 System Generalised Coordinates 

There are two loops in the system, with connectivity 1-2-3 and 1-2-4. Torsion 
spring is not considered as part of loop. Along the loop 1, link 1 is connected 
to the ground (base). Only one orientation is sufficient to define link 1 with 
respect to the base. Hence, fa, i.e. the relative angle between the link 1 and 
base, in the reference system of link zero (ground) is g.c. associated with 
it. Since all joints are revolute joints, g.c. associated with other links are as 
follows : 

* link 2 - fai, relative angle between the links 1 and 2 

* link 3 - fa, relative angle between the links 2 and 3 
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Figure 4.1: Mechanical system with three deformable links 

* link 4 - (j)t, relative angle between the links 2 and 4 
For this system, g.c. vector is 

{X} = [01,02,03»04] T 

The number of g.c. for a link is decided by the type of joint with which the 
link is connected with the previous link. The g.c. for the links connected by 
various joints (with their previous links), are as given below . 

• Revolute joint - 1 (orientation with respect to previous link), 

• Prismatic joint - 1 (displacement with respect to previous link), 

• Spherical joint - 3 (3 Euler angles with respect to previous link), 

• Cylindrical joint - 2 (orientation and displacement with respect to 
previous link) , 

• Screw joint - 1 (orientation with respect to previous link) 
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4.2 Assembly Procedure 


As discussed earlier, providing the correct initial relative coordinates, for 
the complex planar and spatial mechanical systems is a tough task for the 
user. To help him through this, the assembly phase is used. Optimisation 
technique is used to correct the configuration. Fig. 4.2 shows initial and 
corrected configuration for a simple system. All the links are treated as rigid 
links, since there is no deformation at this stage of analysis. The entire 
procedure is explained in this section. 




Corrected Configuration 


Figure 4.2: Correction of user input using assembly 


4.2.1 Joint modeling 

First, from the data provided by the user, various joints in the loop are 
modeled using transformation matrices (refer section 3.2). These individual 
matrices are then concatenated for each loop, so that global coordinates of 
ground connections can be calculated. Using the loop closure equations, the 
coordinates of joints D and E (Fig. 4.1), in the global reference frame are 
calculated as, 

} = T( 0, 00 * T(L U 0a) * T(L 2 , 0 3 ) * {L 3 } (4.1) 
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and 


w = 


E x 

Ey 

E z 

1 


= T( 0, fa) * T(L U fa) * T(L 2 , fa) * {fa} 


These coordinates are needed in formulating constraint equations. 


(4.2) 


4.2.2 Constraint equations 

Number of constraint equations depend on the type of joint at the end of the 
loop. Suppose in the planar system revolute joint is at the end. It creates 
two constraints viz. the global x and y coordinates of the end link at that 
joint, are same as the global x and y coordinates of the ground connection. 

A) Constraint equations for planar systems 

Fig. 4.3 shows the system assembled from approximate data provided by the 
user. Joints D and E are not coincident with the ground connections D’ and 
E'. 


• Revolute joint at the end of loop 

For both the loops, the joints at the ends are revolute. Hence, the constraints 
are that the points D and D' in loop 1, and E and E' in loop 2 must coincide. 
In other words, x and y coordinates of these points must be same. Thus, we 
have, 

{D} = {D'} and {E} = {E'} 

D and E are calculated in previous subsection. D and E are accurate 
coordinates of ground connection. Now, we get two constraint equations per 
loop, which are 

h = D x - D' x = 0, 

; f 2 = D y -D' y = 0, 

,f 3 =E x -E' x = 0, 
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Fi§ur6 4.3i Rcvolutc joint fit the end of a loop 

and 

f< = E y -E' y = 0; 
or 

j £ j = T( 0, *) * r(L„ « . T(U, h) * W - { d; } (4 ' 3) 

and 

| £ |=T(o,*)«T(ti,*)*r(i 1 ,*)*{M-{ £- } ( 44 > 

In tins way, four constraint equations are obtained for this planar system (2 
equations per loop). 
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Figure 4.4: Prismatic joint at the end of a loop 
• Prismatic joint at the end of loop 

The two constraints due to prismatic joint at the end of loop, are that the 
slider moves along a fixed straight line and its orientation with respect to Z 
axis of global reference system is fixed. That is, we have 

y — mx — c = 0 

and 

cj)- (f) c = 0 

cj) c is orientation of slider axis with respect to X axis. Thus, in the Fig. 4.4 
the two constraints obtained for the loop are 

C y - mC x - c = 0 (4.5) 
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and 

4>l + 02 + <^3 ~ <j>c = 0 (4-6) 

B) Constraint equations for spatial systems 

In spatial systems, the number of constraints per loop are three or more 
depending upon the type of joint at the end of loop. Getting the governing 
constraint equations is a bit complicated than planar systems. 

• Spherical joint at the end of loop 

For a spherical joint at the end, we get three constraints : the global X, Y 
and Z coordinates of the end link must be same as the global X, Y and Z 
coordinates of the corresponding ground connection. In Fig. 4.5, position 
of D is calculated using the loop closure equations. Since the user provides 
correct coordinates of D, the three constraints are obtained as 



Figure 4.5: Spherical joint at the end of loop 
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• Revolute joint at the end of loop 

It has one degree of freedom, the orientation about the rotational axis. So, 
it oilers live constraints, three position constraints arc that coordinates of 
points D and D should match as shown in Fig. 4.5. They are obtained sim- 
ilar to that ot the spherical joint. The remaining two orientation constraints 
are obtained as follows. 



Figure 4.6: Revolute joint at the end of loop 

I?.*, is the matrix, we get by traversing the loop from global frame to k - th 
frame. k R m d is the description of end frame with respect to the k - th frame, 
and R. md is the description of the end frame with respect to the global frame. 
So, we have, 


Rk X fc Rend. — Rend (4-8) 

Initially, R md is fixed with its z axis along revolute axis and x, y axis in a 
plane perpendicular to it (forming a right-handed triad). During analysis as 
relative coordinates changes, R md also changes. But, orientation of revolute 
axis in global frame must be same. So, we get two orientation constraint 
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equations by comparing third column (for z axis) of initial R end and updated 
R end during the analysis. 

• Prismatic joint at the end of loop 

A prismatic joint has only one degree of freedom, the translation along the 
given axis (direction), and it offers five constraints. The position constraints 
are obtained as follows. Let l, m, n, be the direction cosines of slider axis. 
Then, its equation is, 

s - s o _ y-y o _ z - zp , . 

I m n \ ’ J 

Now, the reference point on slider should satisfy Eq. 4.9, which gives two 
constraints. Orientation constraints are obtained similar to revolute joint. 
But, since orientation of slider axis is fixed in global frame, we get three 
constraints by comparing initial R end and updated orientation matrix R end . 

4.2.3 Minimisation of objective function 

This is the final step in assembly phase. The constraint equations obtained 
are a set of nonlinear equations. The solution of these equations will give 
us the correct initial configuration. Either we can directly solve the set of 
simultaneous equations, 


K J Tl j 

or, we can pose the problem in such a way that, the sum of squares of all the 
constraint equations is zero, i.e. 

/l 2 + /2 +/f + "- + /» = 0 ( 4 - 10 ) 

Eq. 4.10 can be solved as an optimisation problem 

u = ft + fi + fs + ••• + /« ( 4 - n ) 
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with U &s objective function to be minimized. Thus, it becomes an un- 
constrained nonlinear multi- variable optimisation problem. It is solved in 
MAT LAB with function fmins . The solution obtained is correct initial con- 
figuration, This is analysed for static equilibrium in next section. 


4.3 Static Equilibrium Analysis 

For static equilibrium analysis, major steps followed are : 

• Calculation of lengths of deformable links 

• Generalised force calculation 

• Objective function formulation 

• Minimisation of objective function 

Above steps are explained in detail in the following subsections. 

4.3.1 Calculation of lengths of deformable links 

In static equilibrium analysis, deformation of compliant members must be 
considered. Here, we need to calculate two link lengths, l 2 in loop 2 and 1$ 
in loop 1 (Fig. 4.7). 

First consider link 2 with length l 2 . Suppose, the point B is at (xi,yi) 
and C is at (x 2 , y 2 ), with respect to global reference frame. Fig. 4.7 shows 
the arrangement. Length l 2 is equal to the distance between the points B 
and C. Thus, 



where ykc' 2 + y 2 is the distance between the joints B and C, approached from 
both ground connections of the loop in which the flexible member is present. 

Here, 


X — X-2 ~ Xi 


and 


y - V 2 - Vi 
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Using the loop closure equation, we get, 

{*; }=r(°,«*{£i} (4.13) 

and 

l S } = { e] } - * r (« * T W«) * M ( 4 - 14 ) 



Figure 4.7: Calculation of length of deformable link 

After calculating (xi,j/i) and (; x 2,2/2)) length I 2 can be calculated. In a similar 
fashion, the link length I3 in loop 1 will be calculated. The calculation of 
(xi,yi) and (^2,1/2) is same, irrespective of the type of joint in the loop. The 
only difference it makes is the concatenation of homogeneous transformation 
matrices changes depending upon the type of joints along the loop. 
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In case a deformable member is at the start of the loop, we have, (xi, y t ) = 
(0, 0). Similarly, if the deformable member is at the end of loop, the coordi- 
nates ( x- 2 ,y- 2 ) are the same as the coordinates of the end ground connection 
of the loop. 


4.3.2 Generalised force 

The system deforms under the action of external loading. These forces in- 
clude external forces, torques and body forces i.c. weights of links. 

Fig. 4.8(a) shows external loading acting on the system. Force F is acting 
at a distance L along link 1 , from ground connection. It makes angles a x 
and a y with respect to the global X and Y axis respectively. Moment Mi 
is acting counterclockwise at the c.g. of link 4. Let W' 1 ,W 2 ,W 3 ,W 4 be the 
weights of links 1 , 2, 3 and 4 respectively. First, location of force in a global 
reference frame is calculated. Let, (x, y) be the location of force with respect 
to global reference frame and L = [l x ,l y ,l z , 1] T . Then, 


| * j=r(0, *,).{£} (4.15) 

For Jacobian calculation, the above equation has to be differentiated with 
respect to g.c. The partial derivatives of Eq. 4.15, with respect to 4>y, <j> 2 , 
<j> 3 , (/h are given by 



r'(o,*)*W> 



= { 0 } 


dx 

8<t>z 

JbL 

d<p3 


= { 0 }, 


dx 
d4> 4 
dxi_ 
d<t> 4 


= { 0 } 


Weights of links are acting at c.g. of respective link. Let (xi,yi), (x 2 ,y 2 ), 
(£ 3 , 7 / 3 ), (x 4 ,y 4 ) be the location of c.g. of every link in a global frame. 
{Ci }, {C 2 }, {C 3 } and {C 4 } are their position vectors, from the previous 


joints. Then we have, 


| '2 | = T(0, <f>i) * {Ci} 
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(a) Initial Position 



Figure 4.8: Forces acting on system 
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{ Z }= r (0.A)*r(L„^)»{c 2 } 

{ 2 } = no, <h) * T(L U *) * T(L 2 J 3 ) , { c 3 } 

{ 2 } = no,<h)*T(L h ^ 2 ) , T(L 2 ,<M » {a,} 

The partial derivatives of the above equations with respect to g.c. are calcu- 
lated similar to the procedure discussed above. 

In case of moments, the angle of link on which the moment is acting, is 
expressed in terms of g.c. Let link 4 make angle </> with respect to global X 
axis. Then 


(j> = $ l + <f>2 + $4 

and 

deb r 

dx ~ t 1,1,0,1 ] 

The actual force acting on the system, is the set of force components acting 
along the X, Y and Z directions. Moments are also acting about the c.g. of 
some of the links. The actual force is denoted by F in . 


Fin = [f'"i cos ot x , Fi sin cr v , 0, -V^, 0, -W 2> 0, -W 3 , 0, -W 4 ,+Mtf 


Updating the F in 

As stated earlier, the positions of the external loads are considered fixed, 
with respect to the global reference frame, during the whole analysis. Hence, 
after each iteration, the locations of these forces changes with respect to the 
links on which they are applied. Due to change in position, Jacobian changes, 
which is calculated as given below. 

Consider the force F whose location in global reference frame is calculated 
as (rr, y). With respect to link 1 it is at L{p x ,p y ). After the first iteration, 
the system deforms and hence the location of the point with respect to the 
link 1 changes to (p' x ,p y )- Fig- 4.8(b) shows configuration after iteration 1. 
Suppose, g.c. for link 1 in this iteration is (j>i- Then 
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This gives, 


r (0.#i) *{£'} = | * j 

{L'} = [T(0,^)]- 1 | * | 

Now, the new Jacobian is calculated with position vector L'. Finally gener- 
alised force, F gen is given as 


P gen — J * Fj n (4.16) 

4.3.3 Objective function 

Objective function is the potential function of the mechanical system. The 
factors contributing to potential function are forces acting on the system and 
springs. If there are n springs, with stiffness hi, undeformed length D i} and 
length in a given iteration as di, then their contribution to potential energy 
function is given by 

* = 5 !>(*({*}) -A) J (« 7 ) 

1 1=1 

Torsional springs are assumed to be undeformed at initially correct configu- 
ration. Their contribution to potential energy function is given by 

1 Tit 

^ = U- 18 ) 

z i=i 

Here, n t is total number of torsional springs, k t j is stiffness of j-th spring and 
(O'j — Oj ) is angular deformation. Finally, potential energy due to generalised 

force is 


P3 = 

Sum of the three expressions above gives the objective function <E>. 


(4.19) 


<I> =p i + P2 + P3 


(4.20) 
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4.3.4 Minimisation of objective function 

A mechanical system also contains rigid links, whose lengths are constant. 
Hence, objective function is subjected to constraints as discussed in section 
4.2.2. It is a constrained optimisation problem. Minimisation of is carried 
out by MATLAB function constr , which gives the final equilibrium configu- 
ration. This equilibrium configuration is then graphically displayed. 
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Chapter 5 

Results and Discussions 

The static equilibrium analysis of a few planar and spatial mechanical sys- 
tems is presented in this chapter. 

Presentation scheme is as follows : 

• System input data and numerical output is described first. 

• Below each system input, mechanical system with external forces is 
shown. 

• In the next figure, result of assembly and static equilibrium is shown. 

• Diagram (a) shows assembly output (thick lines) and diagram ( b ) shows 
static equilibrium configuration in thick lines. 

• In dia. (a) incorrect (user-supplied) configuration is superimposed in 
thin lines, and in dia. (6) initial correct configuration is shown in thin 
lines for easy comparison. 
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Example 1 : 


* Mechanical system constitutes of 4 rigid and 2 deformable members. 

* It is a statically indeterminate system and is shown in Fig. 5.1. 

* Link Lengths : h = 10, k = 15, h = 22.8024 

* Stiffness : — GO N/crn k$ = 200 N-cm/rad 

* Ground connection coordinates : A(22. 8024, 0.0) 

* Force : 200 N is acting at revolute joint Ji, Orientation : 90° 

* Moment : 100 N-cm at link 1 

* User input : [40, -GO, 30] 

* Initial correct configuration : [30.0002, -49.4714, 19.4705] 

* Equilibrium Configuration : [42.4762, -68.8051, 26.4288] 

* Final spring length l :i = 20.8195 cm 

* Torsion spring deformation : 0.2177 rad 



Figure 5.1: Example 1 


Mechanical system 
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Figure 5.2: Example 1 : Results 
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Example 2 : 


* System consists of 6 rigid links and one deformable link. 

* It is a statically determinate system and is shown in Fig. 5.3. 

* Link Lengths : l\ = 10, h — h = 10) U = 10 

* Stiffness of torsion spring : 2000 N-cm/rad 

* Ground connection coordinates : A(4.8102, 3.4763), B(23.5364,4.1678) 

* Force : 100 N is acting at slider, Orientation : 0° 


* User input : [100, -100, -120, -30, 50] 

* Initial correct configuration : [114.99, -104.99, -135, -40, 29.99] 

* Equilibrium Configuration : [98.44, -103.54, -144.56, -11.79, 15.08] 

* Torsion spring deflection = 0.2889 rad 



Figure 5.3: Example 2 


Mechanical system 
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(a) Initial Configuration (Ex. 2) 



(1>) Equilibrium Conllguratlon (Ex. 2) 


Figure 5.4: Example 2 : Result 
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• Here Example 2 is analysed for three different forces. 

• Initial and equilibrium configuration is superimposed in Fig. 5.5 


-60 N 



Figure 5.5; Example 2 : Configurations for varying load 

• It is observed that increasing force magnitude increases deformation of 
compliant members. 

• Also increasing spring stiffness decreases its deformation (result not 

shown here). 
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Example 3 : 


* statically determinate system having 4 rigid and 3 deformable members. 

* Link Lengths : l\ = 10, k = 15, k = 9.0, k = 12-5 cm 

* Stillness : k\ « h = 00 N/cm h = 400 N-cm/rad 

* Ground connection coordinates ; A(20. 3541, -0.85251), B(33. 6717, 14.7744) 

* Force ; 600 N at 4 cm from revolute joint J u Orientation : 25° 

= 1 = Moment': 3000 N-cm at link 1 

* Initially slider is at 6 cm along link 2. 

* User input : [35, -10, -45, 30] 

* Initial correct configuration : [30, -20, -60, 25] 

, Equilibrium Configuration ; [40.0361. -27.6175, -90. 12.3189, 10.83] 

* Final spring lengths : I 3 = 9.8454 and k = 12.4646 cm 

* Torsion spring deformation : 0.17518 rad 



Figure 5.6: Example 3 : Mechanical system 
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Example 4 : 


* Mechanical system having 3 rigid and 3 deformable members. 

* It is a statically indeterminate system and is shown in Fig. 5.8. 

* Link Lengths : l x = l 2 = k - U = 10 

* Stiffness : k 3 = k A = 70 N/cm k 6 = 200 N-cm/rad 

* Ground connection coordinates : A(23. 3195,16. 2484), B(19.6592, 2.5881) 

* Force : 150 N at 5 cm from revolute joint Ji, Orientation : —60° 

* Moment : 100 N-cm at link 2 

* User input : [50, 25, -30, -70] 

* Initial correct configuration : [60, 15, -45, -75] 

* Equilibrium Configuration : [47.26, 7.35, -22.95, -84.42] 

* Final spring lengths : l 3 = 9.3595 and Z 4 = 9.7091 cm 

* Torsion spring deformation : 0.22238 



Figure 5.8: Example 4 : Mechanical system 
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(a) Initial Configuration (Ex. 4) 



Figure 5.9: Example 4 : Results 
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Example 5 : 


* Mechanical system having 3 rigid and 3 deformable members. 

* It is a statically indeterminate system and is shown in Fig. 5.10. 

* Link Lengths : h = 10, l 2 = 10, l 3 = 10, PQ = 15, PR = 12 

* Angle RPQ = 60° 

* Stiffness : k 2 = k 3 = 60 N/cm k 5 = 2000 N-cm/rad 

* Ground connection coordinates : A(25.82, 20.58), B (28.32, -1.43) 

* Force : 100 N at 5 cm from A along link 1, Orientation : —60° 

* Moment : 100 N-cm at link 4 

* User input : [20, -40, -30, -50] 

* Initial correct configuration : [30, -44.99, -40, -30] 

* Equilibrium Configuration : [23.1, -46.84, -45.34, -18.98] 

* Final spring lengths : l 3 = 10.56 and I 4 = 9.47 cm 

* Torsion spring deformation : 0.12044 rad 


A 



Figure 5.10: Example 5 : Mechanical system 
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Example 6 : 


* Spatial Mechanical system having 3 rigid and 2 deformable members. 

* It is a statically indeterminate system and is shown in Fig. 5.12. 

* Link Lengths : li = 15, l 2 = 10, U = 12.5 

* Stiffness : = 60 N/crri k 4 = 2000 N-cm/rad 

* Ground connection coordinates : A(30.0, 0.0, 6.0) 

* Direction cosines for axes of revolute joint : 

(0.1227, 0.1227, 0.9848) and R 2 (-0.2645, 0.2438, 0.933) 

* Moment : 90 N-cm at link 1 

* User input : [30, -20, -30, 20] 

* Initial correct configuration : [32.0637, -28.6109, -56.853, 15.7712] 

* Equilibrium Configuration : [34.642, -36.1767, -50.9071, 14.5569] 



Figure 5.12: Example 6 : Mechanical system 
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Example 7 : 


* Spatial Mechanical system having 3 rigid and 3 deformable members. 

* It is a statically indeterminate system and is shown in Fig. 5.14. 

* Link Lengths : li = 10, I 3 — 15, U — 12.5 

* Stiffness : /c 3 = Aq = 60 N/cm k 5 = 1000 N-cm/rad 

* A(15.69, -1.841, 8.604), B(20.248, 1.724, 3.066) 

* Direction; cosines for axis of revolute joint : ( 0 . 1 , 0 . 1 , 0.98994) 

* Force : (135.95, -63.4, 0 ) at 5 cm from A along link 1 

* Moment : 100 N-cm at link 1 

* User input : [25, 5, -45, -35, -40, -25] 

* Initial correct configuration : [29.999, 5, -9.999, -40, -44.999, -19.999] 

* Equilibrium Configuration : [24.4, 4.4788, -41.14, -39.06, -36.14, -19.79] 

* Final spring lengths : Z 3 = 15, h — 11-8765 cm 



Figure 5.14: Example 7 : Mechanical system 
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Bor the mechanical system in Example (7), forces of different magnitude is 
applied. In this case coordinates (X,Y,Z), of endpoint of the link 1, and 
distance of slider along link 1 from ground (L), defines the configuration of 
system completely. They are plotted against force magnitude. 

Stiffness : k 3 = k A = 60 N/cm k 5 = 600 N-cm/rad 



Figure 5.16: Graphs for Example 7 
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Chapter 6 
Conclusion 


6.1 Summary 

In the present work, the static equilibrium analysis of compliant mechanical 
systems is performed. The systems considered are user defined planar /spatial 
systems. Both statically determinate as well as statically indeterminate sys- 
tems are analysed. 

Generalised coordinates used are relative coordinates. The system is mod- 
eled using loop closure equations and joint transformation matrices. In case 
of complicated planar and spatial systems, providing accurate values of rel- 
ative coordinates is a difficult task for the user. Hence, in the first phase of 
analysis (assembly), correct initial configuration of the system is obtained. 
Matlab routine fmins is used for minimisation of objective function. 

In static equilibrium analysis, the compliant members are modeled as 
linear and torsion springs. The external loading is converted into the gen- 
eralised force and potential energy function of the system is formulated. It 
is subjected to constraints that the lengths of rigid links remain constant 
throughout the analysis. This constrained optimisation problem is solved by 
constr function in Matlab optimisation toolbox. Finally, equilibrium config- 
uration in graphical form is displayed. Some results are generated for both 
planar and spatial cases, with various kinematic joints. The results are found 
to be satisfactory in the sense that they do give consistent solutions. 
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Thus, from the results obtained we can conclude that, the methodology 
suggested for static equilibrium analysis is general and computationally effi- 
cient. Main reason for this is the choice of relative coordinates, which leads 
to lesser number of constraint equations compared to Cartesian coordinates. 

6.2 Future Scope 

Now, the analysis is fairly general. It can analyse planar /spatial and stat- 
ically determinate as well as statically indeterminate systems. Hence the 
next step in this research will be to consider combined loading of members 
subjected to torsion, moment and shear along with axial load. Detecting 
singular situations can be another area of research. Also systems analysed in 
present work are idealised. So, considering friction effects at joints will lead 
to more practical approach. 
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Appendix A 


The user input for a sample problem is given below. The sample problem is 
as shown in the figure. The information asked by Software is given in italics. 
Input provided by the user is given in bold type. Explanation is given in 
sans serif type style. 

Enter total no. of links including base and torsional springs : 7 

Is it Planer(l) or Spatial(2) System : 1 

Enter total no. of ground connections : 2 

Enter max. no. of links in loop considering all loops : 4 

Enter loop connectivity matrix : 

3 12 4 10 

4 1 2 5 3 1 

First column contains no. of links in a loop. Subsequent columns contains seri- 
ally link numbers in that loop followed by joint at the ground connection. If any 
column is remaining it will contain zero as entry. 

Joint terminology for 2D : 

Revolute Joint - 1 , Prismatic Joint - 2 
Joint terminology for 3D : 

Spherical Joint - 1 , Prismatic Joint - 2 

Revolute Joint - 3 , Cylindrical Joint - 4 , Screw Joint - 5 

Enter link input matrix : 

1 10.0 0 30 0 

1 15.0 0 -20 0 

2 09.0 0 -60 0 
2 12.5 0 25 0 

5 06.0 0 0 0 
4 00.0 0 0 0 

1st Column : Link Type (1 - Straight rigid link, 2 - Linear spring, 
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3 - Ternary link, 4 - Torsion spring, 5 - Slider) 



Mechanical System 


2nd Column : Link length 

3rd Column : Flag 1 or 2 for ternary link 

4th Column : Relative angle with respect to previous link 

5th Column : Angle between edges of ternary link 

Note : Ternary link is represented by two Straight links with included angle be- 
tween them as constant. 

Enter torsion connectivity : 1 2 

Number of two links which are connected by torsion spring. 

Enter Ground connection coordinates : 

33.671 14.774 
20.354 - 0.85251 

1st Row : X,Y coordinates of 1st ground connection 
2nd Row : X,Y coordinates of 2nd ground connection 

Enter no. of forces and no. of moments acting : 1 1 
Enter input force data : 

2 0 600 25 4 
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1st Column,: Link no., 2nd Column : 0 if not acting on ternary link 

3rd Column : Force Magnitude, 4th Column : Angle w.r.t. positive X axis 

5th Column : Distance from previous link 

Enter input moment data : 1 -300 

1st Column : Link no., 2nd Column : Moment magnitude 

Enter Stiffness of linear springs : 60 60 
Enter Stiffness of torsion spring : 400 



Appendix B 


Functions used in the code for Planar Analysis : 

• sp _len : Function to calculate linear spring lengths. 

• constraint: It models all the joints. Output of this function is ground 
connection coordinates. 

• force_co : Coordinates of point of application of force are calculated 

here. 

• sp_pe : It calculates P.E. due to compliant members. 

• g_force : Here Jacobian, generalised force and P.E. due to it are cal- 
culated. Jacobian is modified in each iteration. 

• All the above functions are called in main program. 

• Command to run main program : x = constr(’ filename’, v, 10) 

• v is initial guess vector of relative coordinates and output vector is x. 

Functions used in the code for Spatial Analysis : 

• joint : Converts g.c. vector to jointwise g.c. 

• sp_con : Calculates link connectivity in loop and finds loop of which 
linear spring is a link. 

• rev_joint : Models revolute joint. 

• sp_len3 : Calculates linear sping lengths. 

• constraint: It models all the joints. Output of this function is ground 
connection coordinates. 
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co_ord_f : Coordinates of point of application of force are calculated 
here. 


• sp_pe3 : It calculates P.E. due to compliant members. 

• g_forcc3 : Here Jacobian, generalised force and P.E. due to it are 
calculated. Jacobian is modified in each iteration. 

• All the above functions are called in main program. 


63 



Bibliography 


[1] D. Livermore, The determination of equilibrium configuration of spring 
restrained mechanisms using (4 x 4) matrix methods, Jr. of Engg. for 
Industry, Feb. 1967, pp. 87-93. 

[2] R. Gupta, Determination of equilibrium configurations of suspension 
systems by a structural deflection method, MS thesis, University of Wis- 
consin, Madison, Wis., 1966. 

[3] J. Garcia de Jalon and R. Aviles, Computer Method for Kinematic 
Analysis of Lower-Pair Mechanisms - 1, Mechanism and Machine The- 
ory, Vol. 16, No. 5, pp. 5f3-556 (1981). 

[4] J. Garcia de Jalon and R. Aviles, Computer Method for Kinematic 
Analysis of Lower-Pair Mechanisms - 2, Mechanism and Machine The- 
ory, Vol. 16, No. 5, pp. 557-567 (1981). 

[5] P. Nikravesh, Computer Aided Analysis of Mechanical Systems, 
Prentice-Hall (1988), 

[6] E. Haug, Computer Aided Kinematics and Dynamics of Mechanical Sys- 
tems, Allyn and Bacon Series in Engineering (1989), 

[7] B. Dasgupta, M. Deshpande and T.Mruthyunjaya, A static equilibrium 
formulation for mechanical systems with compliant members, Mechanis- 
m and Machine Theory (communicated). 

[8] U. Joshi, Static Equilibrium Analysis of Mechanical Systems with Com- 
pliant Members, M. Tech. Thesis, Department of Mechanical Engineer- 
ing, Indian Institute of Tecgnology Kanpur, India (2000). 


64 



[9] R. Aviles, G. Ajuria, V. Gornez-Garay and S. Navalpotro, Compari- 
son Among Nonlinear Optimisation Methods for the Static Equilibrium 
Analysis of Multibody Systems with Rigid and Elastic Elements, Mech- 
anism and Machine Theory 35(2000), pp. 1151-1168. 

[10] D. Hearn, M. Baker, Computer Graphics, Prentice-Hall, (1999). 

[11] J. Craig, Introduction to Robotics, Addison- Wesley, (1986). 


65 



